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We study novel effects associated with electromagnetic wave propagation 

in a Robertson- Walker universe and the Schwarzschild spacetime with a small 
amount of metric stochasticity. We find that localization of electromagnetic 
waves occurs in a Robertson- Walker universe with time-independent metric 
stochasticity, while time-dependent metric stochasticity induces exponential 
instability in the particle production rate. For the Schwarzschild metric, time- 
independent randomness can decrease the total luminosity of Hawking radia- 
tion due to multiple scattering of waves outside the black hole and gives rise to 
event horizon fluctuations and thus fluctuations in the Hawking temperature. 
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I. INTRODUCTION 



Wave propagation and localization in a random media has been studied extensively 
for the last two decades [|2|-|J. It is known that on the mesoscopic scale, classical wave 
propagation in a random media can be treated in a similar way as electron transport in a 
random potential [|]]. One expects to see the diffusive, localization properties of classical 
waves similar to electrons moving in the presence of impurities HH- Wave propagation in 
curved spacetime is an important topic both in general relativity and in semiclassical 
gravity theory ||. Classical scalar, electromagnetic and gravitational waves in Friedmann- 
Robert son- Walker (FRW) Universe probe into the state of the universe and manifest in basic 
cosmological processes such as structure and defect formation, while that in the Schwarzschild 
and Kerr spacetimes depict high energy astrophysical processes in black holes. The second- 
quantized version in terms of quantum fields gives rise to cosmological particle creation 
and Hawking radiation [|l(J which are important processes in the early universe and 
black holes collapse. Recent progress in studying Planck energy processes, especially the 
backreaction effect of quantum fields in curved spacetime, underscores the importance of 
including fluctuations in particle creation [[[]],|T2[| an d the associated energy momentum tensor 
of quantum fields |T^-[r5|, and fluctuations and dissipation in the dynamics of spacetime 

EMI- 

The program devoted to quantum matter field and classical background spacetime with 
metric fluctuations mentioned above is rather involved, because it requires the calculation 



of four point functions |20[ and demands a self-consistent solution |2T],^2]. In this paper, 
as a useful parallel, we attempt to address an easier problem, that of wave propagation 
in a stochastic spacetime. It is designed to highlight the effect of fluctuations in a back- 
ground metric, while not demanding an explanation of their source, or their mutual influ- 
ence. Stochastic components in the metric can be induced by primordial gravitational waves, 
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topological defects in the sub-Planckian scale, or intrisic metric fluctuations of background 
spacetimes at the Planck scale. Their detection and analysis can provide valuable informa- 
tion about the state of the early universe and black holes. As distinct from the self-consistent 
treatment which is necessary for Planck scale processes, wave propagation in curved space- 
times with metric stochasticity is a test-field treatment. It is nonetheless still a useful probe 
for fluctuations in sub-Planckian processes (such as GUT scale phase transition), which could 
have left important imprints on the observable universe |23fl . 

In this paper, we wish to study novel effects associated with electromagnetic waves prop- 
agation in the Friedmann-Robertson- Walker (FRW) universes and the Schwarzschild space- 
time with a small amount of metric stochasticity. Here we employ a useful observation to link 
up with the more familiar subject of wave propagation in random media studied extensively 
in condensed matter and mesoscopic physics |||| . We first show the formal equivalence of 
the wave equations in curved spacetimes with wave propagation in a media in flat space 
and identify how the metric components appear in the permittivity function (or refractive 



index) of the media |24|. Then we introduce metric fluctuations as a stochastic component 
in the permittivity function and study wave propagation in a curved spacetime with metric 
stochasticity as if it were in a random media. In a spherically symmetric spacetime the wave 
equation for the radial part can be written in the form of a Schrodinger equation in one 
dimension. The effect of the curvature of spacetime appears in the potential term in the 
equation. Once the wave equation is reduced to the parametric form with a stochastic com- 
ponent, the familiar methods used in quantum field theory in curved spacetime and insights 
accumulated in mesoscopic physics can work in each other's advantage. We analyse in de- 
tail the cases with time-independent and time-dependent metric stochasticity and find that 
localization of electromagnetic waves occurs in a metric with time-independent fluctuations. 
In cases where there are time- dependent randomness in the metric, exponential instability 
in the particle production rate occurs. These are new effects due exclusively to the presence 



of metric fluctuations in the background spacetimes. For the Schwarzschild metric, time- 
independent randomness outside the horizon will decrease the total luminosity of Hawking 
radiation due to multiple scattering of waves, if the randomness reaches the horizon, it 
contributes to the fluctuations of Hawking temperature. Time- dependent stochasticity is a 
more complicated matter which requires a self-consistent analysis of the interaction between 



spacetime and the waves or fields as encoded in the fluctuation-dissipation relations |[21| , [22 
This is to be investigated later as a part of the stochastic backreaction problem mentioned 
above. 

The paper is organized as follows: In Sec. 2, we show how an electromagnetic wave 
propagating in curved spacetime can be related to that in flat space but with a refractive 
index depending on the metric components ||24|| . This section is meant to be a short cut for 
readers not too familar with curved spacetime physics to see the correspondence with wave 
phenonema. (Readers familiar with it can skip to Sec. 3). In Section 3, using the methods 
developed in Sec. 2 we study wave propagation in curved spacetime where the metric has 
a stochastic component. We use the Friedmann- Robertson- Walker universe as an example 
to show that if the stochastic component of the metric is independent of time, and for a 
sufficiently smooth randomness the Maxwell equation has the same form as a conformally 
coupled scalar wave equation. If the randomness appears only in the radial direction, the 
radial wave equation has the same form as an one dimensional Schrodinger equation in a 
random potential. Because the electromagnetic wave equation is conformally invariant, it 
bears the same form as in flat space. With a stochastic metric component we see that wave 
localizes in space. In Sec. 4, we study the case of time-dependent (but space-independent) 
stochasticity in the metric, and show that parametric amplification takes place giving rise 
to cosmological particle creation. However, because of the metric fluctuations the rate of 
this amplification increases exponentially in time. The fluctuations in the particle creation 
rate is also discussed. In Sec. 5, we study wave propagation in a Schwarzschild spacetime 



with time-independent fluctuations. We show that this can decrease the total luminosity of 
Hawking radiation due to multiple scattering of waves outside the black hole. Fluctuations 
which involve event horizon will give rise to fluctuations in the Hawking temperature. In 
Sec. 6 we summarize our findings and end with a short discussion. 

II. CLASSICAL ELECTROMAGNETIC WAVES IN CURVED SPACETIME 

Maxwell's equations for an electromagnetic field tensor F^ u in a gravitational field with 
metric g^ u are given by 

F^ = F^, a + F u ^ + F^, u = {111) 

where semicolons denote covariant derivatives with respect to g^ v . These equations can 
be cast in a form for waves propagating in a permeable media in flat space ||24|| . The 
correspondence between E, H (the electric field and magnetic induction) and D, B (the 
electric displacement and magnetic field) respectively is given by 

Di = €ikEk — {G x H)i 

Bi = fi lk H k - (G x E)i {11.2) 

where the dielectric and permitivity functions are 

/ \i/2 9 lk 

tik = Hik = -{-g) 1 — 

Oi 

Gi = - 9 —. (II.3) 

#00 

A. Friedmann-Robertson- Walker Universe 

Let us now consider the Friedmann-Robertson- Walker (FRW) spacetimes with line ele- 
ments 
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ds 2 = -dt 2 + a(t) 2 R 2 [d X 2 + s 2 ( X )(d6 2 + sin 2 6d<j) 2 )} 



(II.4) 



where s(x) = sinx, x, sinhx correspond to closed, flat, and open cases respectively. Here 
a(t) is the scale factor and R is the radius at time t where a(t ) = 1. 
Using Cartesian coodinates, 

dt = a(t)dr) 

x 1 = 2Rt(x/2) sin 9 cos 
x 2 = 2Rt(x/2) sin 9 sin 

= 2Rt( X /2)cos6 (II.5) 

where rj = J dt/a is the conformal time, and t(x/2) = tan(x/2), x/2, tanh(%/2) correspond 
to the closed, flat, and open cases respectively. We can then write the FRW line element in 
the form 



- ds 2 = a 2 [~dr] 2 + f 2 (p)(8 ij dx i dx j )} 



(II.6) 



where 



p = E(-T] 1/2 

i 

= 2Rt( X /2) 



(11.7) 



and 
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\ l-p 2 /4R 2 / 

where the column elements correspond to closed, flat, and open cases respectively. 
In this isotropic form of the metric, the Maxwell equations are given by 
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V x E = -^(f(p)H) 

V • (fE) = 

V • (fH) = (II.8) 
where one can identify the dielectric and permittivity functions as 

Qfc = Pik = f{p)$ik 

d = 0. (II.9) 

Assuming a harmonic (conformal) time dependence of the solution e~ luJ ' n , we can write 
( gig ) as 

VxE = ioof(p)H 
VxJJ = -iujf(p)E 

V • (fE) = 

V •(/#) = (). (11.10) 
The symmetry in the relation 10| ) allows us to write this in the more compact form 

V x F = uf(p)F 
V-(/F) = (11.11) 

where F = E + iH. 

For the FRW metric with spherical symmetry the solutions can be expressed in terms 
of the vector spherical harmonics Yi m (9,<j)) obtained by operating on the scalar spherical 
harmonics Yi m (6, <fi) with the invariant operator — ir ■ V as 



E = Et^f m V x gf(p)Y lm + AZgf{p)Y l 



, uf 

Lm J 



lm\ 



H = Y\ A ?m9?{p)Yim ~ —Ay x gr{p)Y lm ] (11.12) 

l,m 
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where gf' M {p) are functions of radial distance only and Af^ M are the coefficients. Super- 
scripts E and M denote the electric and magnetic multipole field components, respectively. 
The field F is given by 

F = EMk(V x gf(p)Y lm + gf(p)Y lm ) 

l,m U J 

+ Afi(gf(p)Y lm + -^V x gf\p)Y lm )}. (11.13) 
One can identify the (I, m) component of a magnetic multipole field as 

F% = g^{p)Y lm + -i-V x gf{p)Y lm . (11.14) 
Writing the magnetic multipole scalar field function 

^ M (x) = Rp{x)g?\p(x)), (H.15) 

we obtain the radial equation for ip M as 

rf V%) , , a rT , \\.i,M 



d X 2 



where n = Rw and 



+ (n'-U(xW V1 (x) = (H.16) 



tf(x) = (n.i7) 



B. Schwarzschild Spacetime 

The line element for the Schwarzschild spacetime is given by 

If we define r = p + M + M 2 /4p, we can write ( [111.35] ) as 

- rf s 2 = -f x (p)dt 2 + ft(p)[dp 2 + pW + p 2 sin 2 #d0 2 ] (11.19) 



where fi(p) — 1 — and f2(p) = f/p. Furthermore, a spatial coordinate transformation 

x 1 = p sin 9 cos <p 
x 2 = psin#sin0 

x 3 = pcos9 (11.20) 

brings ( |II.19| ) into a Cartesian form 

- ds 2 = -fi(p)dt 2 + fl{p)[{dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 }. (11.21) 

Considering this spacetime acting as a medium, the corresponding dielectric and permittivity 
functions are expressed in the same form as in ( |II.9| ) with 



In such a medium, the radial equation for the magnetic multipole field is [25 



d 2 jj M (r*) 
dr* 2 



+ [u 2 - U(r*)]i) M (r*) = (11.23) 



where r* = r + 2Mlog(r/2M — 1) is the Regge- Wheeler (turtoise) coordinate, and the 
potential term is given by 

f(»H^±->-^]. (H.24) 

r z {r*) r{r*) 



III. WAVE PROPAGATION IN CURVED SPACETIME WITH METRIC 

STOCHASTICITY 

Consider the metric as containing an averaged (deterministic) part gffj and a small fluc- 
tuating (stochastic) part 



9,u = 9 { S + K u (III.l) 
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then from ( |il.3| ) 



'_ o (0)U/2 t 

- 9 > .g(0)ik + 



-g(0)V/ 2 



(0) 
#00 



(0) 
#00 



4? + Se ik 



(III.2) 



where we have used the synchronous h 00 = gauge |26| . The traceless h = Z)f=i h\ = gauge 
further simplifies the expression and we obtain the correlation function of the fluctuating 



part of the refractive index in terms of the metric fluctuations as 



(Se ij (x)Se k i(x')) 



(-^(0)( :E )) 1 /2(_^(o)( x /))i/2 



-(h l \x)h kl (x')) 



For flat spacetime, it has the simple form 



(III.3) 



(Se^Seuix 1 )) = (h^(x)h kl (x')). 



(III.4) 



We see metric fluctuations can thus be represented as fluctuations in the optical index. 



A. Time Independent Metric Stochasticity 

In the rest of this section, we consider cases in which the permitivity function contains 
only spatial disorder. 

If we assume the background metric as well as the stochasticity are both isotropic, we 
can write the corresponding refractive index as 

Gi = (III.5) 

where f(x) is a random variable. 

Assuming a harmonic time dependence of the field with frequency lu, Eq. ( |11.11|) with 
source term J becomes 
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{L — tof(x))P = J 



(111.6) 



where L = Vx . Define Green's function g(x,x f ) as 

(L - uf(x))g(x, x') = 8(x - x'). (III. 7) 

Then the disorder-averaged Green's function is 

G(x,x') = (g(x,x')) 

= ((L-uf(x))- 1 }. (111.8) 

If we define the mean field Green's function (x, x') as 

G {0 \x,x')^(L-oo{f(x)))- 1 (HI.9) 
we obtain the Lippman-Schwinger equation 

G-\x, x') = G {0) -\x, x') + wS(i, x') (111.10) 
where E(x, x') is the self-energy. Also averaging (|III.6|) , we have 

(LF)-u(f(x)F) = J. (111.11) 
Define f = fF = D + iB, then from ( pTTTI) 

u{P) = {IP) - J 

= J dx'(G^- l (x,x') +ujZ(x,x') -L)(P) (111.12) 

(P(x)) = f dx'((f(x))5(x-x') -E(x,x'))(F(x')). (111.13) 



We see that the complex electromagnetic displacement vector has a nonlocal dependence 
on the media due to the self-energy term. We can thus develop a formalism making use 
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of techniques from quantum field theory to study the effect of dissipation, higher order 
correlation, screening, etc. (see e.g., p7||). This is a subject for future investigation. 

Further simplification arises in cases when the polarization effects are negligible. From 
( |11.10| ), we obtain 



V 2 E + u 2 f 2 (x)E + V log /(f) x Vx E + V((E- V) log /(f)) = 0. (111.14) 

Under the conditions that the inhomogeneity and disorder in / are smooth, namely, ujR » 1 
and u\ » 1, where A is the characteristic length scale of the disorder, the third and the 
fourth term representing polarization effects are negligible. (|III.14|) then reduces to 



V 2 E + u 2 f 2 {x)E = 0. (111.15) 

Then each component satisfies the scalar wave equation 

V 2 (j) + uj 2 f 2 (x)(f) = 0. (111.16) 

Now consider the conformal type of stochasticity only in the spatial part of the metric of the 
form [] 

/(f) = f(p)e^\ (111.17) 

where cx(f) is a stochastic function and g is a small parameter. Expanding it with respect 
to g as 

/(f) = /(p)(l+^(f)), (111.18) 

IBT) becomes 



The conformal invariance of Maxwell's equations implies that the conformal fluctuations in the 
spacetime metric can be factorized away, and hence they do not give rise to particle creation [28]. 
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V 2 + oo 2 f 2 (p)(j) + 2gu 2 f\p)a{x)(j) = (111.19) 

In a flat FRW spacetime f(p) = 1, or, for p « R in closed and open cases, f(p) ~ 1. In 
this region, all three cases are described by the equation in flat spacetime 

V 2 + uj 2 ct) + 2guj 2 a(j) = 0. (111.20) 

The corresponding equation for the Green function is 

(V 2 + uj 2 + 2gu 2 a)G(x, x', uj) = 5(x - x'). (111.21) 

The disorder-averaged Green function (G(k, k',u)) = G(k,u)5kk' is known to have a form 

G(k,u)= I (111.22) 
cu^ — k 2 + 2 j 

in the leading order in (ul)^ 1 , where I is an elastic mean free path. For a Gaussian white 
noise (local) random variable ct(x) , / = 2n/gDuj' 1 where D is the strength of disorder defined 
by (a(x)a(x')) = D5(x — x). 
We also define the intensity 

I(x) = (<f>*(x)<f>(x')) (111.23) 
and the intensity-intensity correlation function 

C(x,x') = (<f)*(x)(f)(x)<j)*(x')<j)(x')}c (111.24) 
where C denotes the cumulant. As shown in [|29||, we can evaluate the disorder-averaged 



intensity as 

«*» = tJ^- (HI.25) 

If we confine the disorder in a slab of thickness L, and cross-sectional area A, then for 
r = \x — x'| < /, 
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C(x,x') = (I(x))(I(x')}(^^-) 2 e~T. 



ur 



(111.26) 



For r > I, diffusion modes interact with each other and give rise to a long range correlation 
with a power law decay f30fl . In the Schrodinger equation analogy, this corresponds to 



conductance fluctuations in metals [311. 



B. Radial Disorder — Localization 

Here we also start from the conformal type of metric stochasticity but assume that the 
stochasticity is only in the radial direction 

f( x ) = l^^e 9a{x) (111.27) 

where g is a small parameter ( Such kind of randomness may arise from averaging out 
more inhomogeneous types of randomness in a small spatial region due to the spherically 
symmetric nature of the underlying curved spacetime). 

If <j{x) varies slowly in space (V « a ), the wave equation in the radial direction ( |II.16| ) 
is replaced by 

^ + [n\l + gaf - = (HI.28) 

where we simply use ip for ip . This has the same form as a scalar wave equation in a media 
with a random refractive index. 

For the closed FRW universe, the angular-momentum induced potential term in (|IH.28| ) 



is a single well potential which becomes infinitely high at the origin and x = 71 ■ Then the 
wave is necessarily localized around the bottom of the well in this coordinate. On the other 
hand, for the flat and open FRW universe cases, the potential term becomes asymptotically 
flat as the radius becomes infinite. In these cases, the stochastic term plays a significant 
role in the transport property of the wave. If we expand the second term in ( |111.28| ) into 
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n 2 + 2gn 2 a around the perturbation parameter g and note the positivity of the original 
second term, we see ( |I11.28| ) describes a Schrddinger equation with a fixed positive energy. 
As such , we can write ( |ill.28| ) as 

" jt + V (XW = W ( m -29) 

where E = n 2 and V(x) = —2gn 2 a + + l)/s 2 (x)}- 

For the flat or open universes, the potential barrier vanishes in the asymptotic limit. 
The property of the eigenfunction in this case is well-studied from the context of electron 
transport j§-[(J . It is known that in one dimension for any value of E , the eigenstate of Eq. 
( |U1.29| ) localizes P],[|]. Namely, all eigenfunctions decay exponentially with rate given by the 



Lyapunov exponent [see ([III.32| ) below]. If the correlation radius \ c of the fluctuating part 
of the potential V(x) is smaller than the "wavelength" f\ (x c << n the potential can be 
considered as the white noise type 

(V(x)V( X ')) = Ag 2 n\a{ X )<x)) = Ag 2 R^D5{ X - x!) (HI.30) 

where we have assumed a delta function type of correlation with the strength D for the 
random potential <t(x). 

The cumulative density of states N(ou) has the following form for the white noise potential 
with correlation given in ([111.30 ) p2| . 



(2R 2 u 2 D) l l z 1 
N ( UJ ) = ~2 A^i-iRuo^D) 2 ^) + Bi 2 ( — (Ruj/2D) 2 / 3 ) (1U-31) 

where Ai and Bi are two independent Airy functions which satisfy the equation y" — xy = 0. 

The asymptotic behavior of the solution shows an exponential growth characterized by 
the Lyapunov exponent defined by 



2 The "energy" E = R 2 oo 2 in ( III.29|) is dimensionless and so is the angular variable x- 

3 Hereafter we absorb g in the definition of D 
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X w = lim XM = OT 1 ( s (HI.32) 

x-oo 2Li oc [(j) 

where Li oc (uj) is the localization length given below and the wave function grows as 

tfj( x ) -> e A " (x)x (111.33) 

as x increases. 

In the case of white noise potential given in ([III.30I ), the Lyapunov exponent has the form 



A. = ^>(w) f ^e-y'/^/^ydy. (111.34) 

The entire curve is plotted in Fig.l. We obtain the asymptotic behavior of the localization 
length 

[ (D/Rijj) 2 ^ for u -> 
L loc {u) ~ (111.35) 

1/D for uj — > oo. 

We see that in the long wavelength limit uj —>■ 0, the wave function delocalizes. This is a direct 
consequence of the multiplicative nature of the stochasticity for classical wave propagation 
in random media |3J . 

IV. WAVE PROPAGATION IN SPACETIMES WITH TIME-DEPENDENT 
METRIC STOCHASTICITY: PARAMETRIC AMPLIFICATION AND PARTICLE 

CREATION 

Due to the conformal invariance of Maxwell's equations, time-dependent fluctuations in 
the scale factor can be transformed away (in a conformally- related coordinate). Assuming 
spatial homogeneity ,the time dependence of the metric is uniform in space, we start from 
the metric with the conformal fluctuations similar to ( |III.17|) , but let the spatial part of the 
metric / acquiring a time dependence. 
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If the background spacetime curvature is negligible, the corresponding dielectric and 
peameability functions take the form 

£%k = IHk = h(rj)S ik 

Gi = (IV. 1) 

where we write the time-dependent factor as h(rj) instead of f(rj) for clarity. Assuming h(rj) 
has the exponential form as in (|III.17| ) 

h( v ) =e u(ri \ (IV.2) 

then Maxwell's equations become 

V x E = -JifhMH) 

drj 

V x H = ^-(h(rj)E) 

drj 

V-E = 

V • H = 0. (IV.3) 

From this we obtain a wave equation for E 

VxVxi= -^LfhM-^hME), (IV.4) 
drj drj 

and a similar one for H. 

Since E has no divergence, for a slowly varying noise such that \h(j])\ << ujh(i]), where 
io is the frequency of E, the above equation takes the simpler form 

S^-ii^ = ° < IV ' 5 > 

In momentum space, it reads 
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Expanding in terms of v gives 

^E k + [k 2 - 2u( V )k 2 ]E k = 0. (IV. 7) 

This equation should be compared to ( |111.29| ) with j = by identifying V = 2vk 2 and 
E = k 2 . 

Following the argument from (|IH.29| ) to (|III.35|) and noting that the time coordinate 



r] plays the role of space coordinate \ here, we see the qualitative features of solutions to 
( [IV. 7|) : they show an exponential parametric amplification |3t| with the rate characterized 
by the Lyapunov exponent A^ introduced in ( |III.32|) 

A, = lim \ h ( V ) = — i- (IV.8) 
v-^oo 2L h (k) 

and E k grows 

E k -> E k ( v = 0)e x ^)v (IV.9) 

as 7] increases. The characteristic time scale for parametric amplification Lh(k) is equivalent 
to the localization length in ([IH.32j ) under the rj — > x correspondence. For a Gaussian white 
noise v with correlation (K^KV)) = D5(i] — ?/), the asymptotic behavior of Lh{k) is given 
by 

I (D/k) 2 / 3 for k -> 
L h (k) ~ (IV.10) 
I 1/D for k -> oo. 

It is well-known that the cosmological particle creation problem in a FRW universe can 
be cast into an one-dimensional wave scattering problem || by reading the time variable r\ 
as a space variable, and adopting the following boundary conditions: 



-ikr) 



E 



ae 



-ikr) 



for 7] — » — oo 



In this analogy, the reflection coefficient (3 (or rather, |/5| 2 )in the wave scattering picture 
gives the parametric amplification factor and, in second quantization, the particle creation 
rate. Since in this picture the transmission coefficient t corresponds to the incoming flux 
which is normalized to one, the effect of localization translates to an exponential increase of 
the particle creation rate (3 (so is a) 

P^ e *h(n)v_ (IV. 11) 

More explicitly, if we compare ( IV.ll ) with the corresponding scattering problem in which 



the wave packet is incident from rj = oo and transit to r\ — — oo 



E 



te lkr] for r] — > — oo 

e -ikri _|_ re ikr) f Qr ^ _^ ^ 

we see that the particle creation rate N = \/3\ = \r/t\ 2 in expanding universe measures 
a dimensionless resistance in the wave transport picture. The transmission coefficient in 
a random potential is a multiplicative quantity which decreases exponentially with length. 
This implies that the logarithm of the resistance is an additive quantity and its average 
in length becomes constant in the thermodynamic limit and hence non-fluctuating. The 
average of the resistance itself (hence the average particle creation rate) is also known to 
show the exponential increase in length ( in conformal time rj ) with the rate characterized 
by Aft. We can obtain the asymptotic behavior of the average of N in the high energy limit 
as follows 

(iV) = <|/5| 2 ) = I(e 2 ^ + 1). (IV.12) 
Hence the fluctuations of (N) becomes 

AN 2 = (N 2 ) - (N) 2 = l e *v/Lh( e 2v/L h _ !) ( IV .13) 
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for 7] » Lh where = l/D is the characteristic particle creation rate which appeared in 
( |1V.8| ). The exponential divergence of AN 2 results in the non-self-averaging nature of the 
observable. 



V. HAWKING RADIATION FROM A FLUCTUATING BLACK HOLE 

In this section, we study the effect of metric stochasticity in the Schwarzschild spacetime 
and discuss how it affects the Hawking radiation JlO| . 



In the presence of metric stochasticity, the Helmholtz equation acquires a fluctuating 
component in the refractive index similar to the Friedmann-Robertson- Walker (FDR) uni- 
verse case. Assuming the same conformal form of stochasticity as in ( |III.17|) , but in a radial 
direction only for simplicity, the Helmholz equation in this case has the following form 



d 2 ^ M (r*) 



+ K + 2a(r*)u 2 - U{r*)}ij M {r*) = (V.1) 



dr* 2 

where a(r*) is a stochastic variable as a function of r*. If we assume the metric stochas- 
ticity is restricted in the finite region — L < r* < 0, the incoming flux from future infinity 
suffers backscattering not only from the potential U(r*) but also from the randomness. The 
absorption probability is thus reduced by a factor of e~ x ^ L where A w is the Lyapunov 
exponent which appeared in ( |1H.32|) with R = 1. 



Hence the total luminosity has the form 

1 OO fOD 

L = —J2(2j + 1) d^T e- x - L /{e^ - 1). (V.2) 
2vr ~ Jo 

The more interesting case is when the fluctuating region contains the event horizon such 



that the stochasticity induces horizon fluctuations p^-|36|. For simplicity, we restrict the 
arguments to the scalar wave case in the rest of the section. The scalar wave in the metric 
( |11.19| ) satisfies the equation 
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This takes on a simple form if we define the coordinate x such that dx = f(p)dp/p 2 

[|^+aA 4 (x)]0(x) = O. (V.4) 

If the metric in ( |11.19| ) has the conformal type of stochasticity in the form fi(p)e 2a ^ and 
f2(p)e u ( p \ where <r(p) and z/(p) are two independent random variables with zero mean, the 
equation corresponding to ( |V.3| ) becomes 

|[pr(p)(l - ^)VV^>-"M^ + „V(p)j 0(p) = 0. (V.5) 

We can absorb er(p) by choosing the coordinate x such that = f (p)e u ^~ u ^ dp / p 2 . Then 
we have 

i + uJ 2 r\x)e Au{p) U{x) = 0. (V.6) 
dx 1 

At the horizon, only incoming wave exists. So the solution of ( |V.6|) will be 

~ exp[-iujr 2 e Mps) x} (V.7) 
where r s = 2M and p s = M/2. This mode function behaves as 

-> 1 - ia;r 2 e Mps) x 

^ x + !^ (v _ 8) 

p 

as x — > 1 or, equivalently, p — > oo. Here we assumed the stochastic variables er(p) and z/(p) 
vanish at infinity 

Far from the horizon, the asymptotic flatness of the spacetime allows us to write ( |V.3| ) 

as 

|[pr(p)^] 2 + ^r 4 (p)|0(p)=O (V.9) 
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with the solution 



up 
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[AJ 1/2 {ujp) + BJ_ 1/2 (up)] 
for up « 1. 



up 



(V.10) 



Comparing (|V.10| ) with (|V.8| ), we obtain the following absorption coefficient for low energy 
scattering of the scalar wave, 

' l + i- 



1 - 



1 1 A 



4cA 2 e Mpa) for uo -> 0. 



From ( |V.11| ), the absorption cross section for this process is 



(v.n; 



71 

r 

n 1 a 



47rr s 2 e 4l/(Ps) for -f 
= 47rr s 2 [l+4z/(p s )+8z/ 2 (p s ) + - 

The average absorption cross section is then 

(a a )=4nr 2 s +32nr 2 s (u 2 (p s )) 
= A H + 8A H (u 2 (p s )) 



(V.12) 



(V.13) 



where Ah = Airr 2 is the area of the event horizon and assume all the higher moments of 
v vanish. It is known that low energy absorption cross section of the scalar wave gives the 



area of the event horizon for a large class of black hole solutions [p8 



Therefore, we expect the second term gives the horizon area fluctuations 



SA H = 8A H (u 2 ( Ps )) 



which induces fluctuations in the Hawking temperature 



(V.14) 



(V.15) 



-H 
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VI. DISCUSSION 



In this paper, we have studied the effect of metric fluctuations on wave propagation in flat 
and curved spacetimes. We saw that for electromagnetic waves the effect of metric stochas- 
ticity is equivalent to that of a random optical index. With this analogy, we can calculate 
the correction to the scattering cross section of the electromagnetic wave in curved spacetime 
with metric fluctuations. We also see the intensity fluctuations due to multiple scattering. 
If we assume classical spacetime is effectively emerging from the underlying quantum fluctu- 
ations by coarse graining, it seems reasonable to assume that the stochastic part reflects the 
symmetry of the background spacetime. When the stochastic part of the metric has only 
radial dependence, the problem is further reduced to an one dimensional transport prob- 
lem in which Anderson localization is manisfest. For time-dependent but space- independent 
stochasticity, particle creation rate shows exponential instability in conformal time. Fur- 
thermore, analogy with mesoscopic transport in one dimension predicts large fluctuations in 
the particle creation rate due to its non-self-averaging nature. For Schwarzschild spacetime, 
time-independent stochasticity induces fluctuations of the event horizon and , correspond- 
ingly, Hawking temperature. 

The restrictions we made on separate space and time dependence of the stochastic metric 
is mainly for technical simplicity. We expect qualitatively similar phenomenon in the more 
general cases, where space- and time -dependent metric fluctuations both exist. Certain 
modifications are expected such as dissipative effects on localization for a time-dependent 
random potential. 

Possible origins of the metric fluctuations have been discussed recently. Primordial 
stochastic gravitational waves is currently under intense investigation p^|39|j40| . Stochas- 
ticity used in such a context is different from that in stochastic semiclasical gravity theory 
|12] , p!6| -[18| . In the former it refers to (quantized) linear perturbation of the metric obeying 
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Einstein's equations whereas in the latter, metric fluctuations is induced by fluctuations in 
quantum matter fields, which obey the Einstein-Langevin equation, as a generalization of 
the semiclassical Einstein equation governing the mean values of metric and matter stress 
energy. The naive replacement of the stress-energy tensor with its expectation value leads to 
some pathological results such as the violation of the weak energy condition of the Einstein 
equation [|nj. Possible cure may be obtained by including the smearing term in the classical 



background spacetime |L5|. Self-consistency between the classical gravitation and the quan- 
tum matter sectors also points to the inevitable dynamic role of metric and field fluctuations 



11]. From the astrophysical viewpoint, cosmic string network may provide the source of the 



stochastic gravitational wave background which lies in the observable frequency range of the 



LIGO detector p3 |. The discussions and results in this work are more of a generic nature 
insensitive to the origin of particular sources. 

Ford argued that the lightcone fluctuations based on linearized quantum gravity induce 



black hole horizon fluctuations [[34]]. He suggests that the horizon fluctuations are small 
enough to validify the semiclassical derivation of Hawking radiation. Casher, et. al. argue 
that the semiclasscal approximation breaks down close to but much larger than the Planck 



scale from the horizon |36| due to the interaction with the atmosphere of the horizon. Sorkin's 
argument is based on Newtonian mechanics and its vadility for general relativistic cases re- 



mains to be shown |35[ . Since we assumed a static metric stochasticity for the Schwarzschild 
spacetime, direct comparizon with their results is not obvious. Nevertheless, we believe our 
arguments have certain degrees of generality, knowing that the randomness can arise from 
the imhomogeneity of a static spacetime. These issues are currently under investigation. 
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